The complex channel networks of bone structure 
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Bone structure in mammals involves a complex network of channels (Havers and Volkmann chan- 
nels) required to nourish the bone marrow cells. This work describes how three-dimensional recon- 
structions of such systems can be obtained and represented in terms of complex networks. Three 
important findings are reported: (i) the fact that the channel branching density resembles a power 
law implies the existence of distribution hubs; (ii) the conditional node degree density indicates a 
clear tendency of connection between nodes with degrees 2 and 4; and (iii) the application of the 
recently introduced concept of hierarchical clustering coefficient allows the identification of typical 
scales of channel redistribution. A series of important biological insights is drawn and discussed. 



Containing two types of channels, namely Havers and 
Volkmann channels, the bones of mammals are pervaded 
by a network of intercommunicating channels required to 
maintain the bone living cells jj] , which is done by blood 
vessels distributed inside those channels. Although play- 
ing a decisive role in animal development, diseases, and 
bone reconstitution, the topology of the bone channel 
network remains a subject of investigation. Because of its 
generality, the concept of complex networks |3, i3] stands 
out as a natural and comprehensive means to represent, 
characterize and model the bone channel system. By as- 
sociating the branching points of channels to the nodes 
of the network, and representing the connections between 
such spots in terms of arcs between the network nodes, 
it is possible to represent the bone channel structure in 
terms of a complex network. Topological measurements 
obtained from such networks — especially the node de- 
gree density 0, the conditional node density [5|, and 
the hierarchical clustering coefficient 'g!, Q — can then 
be used to identify important underlying features of the 
analyzed channel system. 

The bone used as a biological model in this study is 
the humerus of the adult cat, collected from necropsy 
cats in the Sector of Pathology at the Veterinary Hospi- 
tal of the Federal University of Uberlandia, Brazil. After 
dissection of the skin and muscles, a ring 0.5 cm wide 
was removed from the central portion of both right and 
left humerus of each animal. Such rings were fixed in 
10% formalin for at least 48 hours. After setting, these 
samples were decalcified in 4% nitric acid per 30 days, 
embedded in paraffin according to a classic histological 
technique, and sectioned in microtome in order to obtain 
200 serial sections with 5 /im thickness. These sections 
were then stained by the Schmorl procedure, and digi- 
tal images were acquired by using an Olympus Triocular 
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BX40 microscope connected to an Oly-200 camera inter- 
faced to a PC-compatible computer through digital plate 
Datatranslation 3135. After the digitalization procedure, 
we extracted the regions of interest through the segmen- 
tation 's'l of both Havers and Volkmann channels. 

The images with final size of approximately 700x480 
pixels were converted into PNG format and used to ob- 
tain a three-dimensional reconstruction of the studied 
structure, which was performed by using the vtkPN- 
GReader class in the Visualization Toolkit for C4-+. 
Figure H^a) shows the three-dimensional reconstruc- 
tion j§, 10] of part of the analyzed bone structure (one 
fourth) . The detection of nodes of the graph and respec- 
tive connections was performed manually, yielding an un- 
weighted, non-oriented complex network, which was rep- 
resented in terms of its iV x adjacency matrix K, such 
that the existence of a connection between any two nodes 
i and j is represented by making K{i,j) = K{j,i) = 1. 
The vertices of "V"-like channels were understood to cor- 
respond to branching spots with missing connections, 
giving rise to nodes with degree 2. Figure ^b) shows 
the geographic complex network (the nodes have posi- 
tions in R^) obtained from the structure in Figure QJa) 
by projecting the node coordinates along the z— axis. 
The degree of a node i is defined as the number of 
edges connected to that node, which can be calculated 
as k{i) = J2q=i ^(*' l)- FigureEIa) shows the dilog plot 
of the node degree density obtained for the bone consid- 
ered in the present study, which is characterized by 988 
nodes and 1120 links. 

The obtained density profile indicates that the bone 
channel network resembles a scale free model , which 
is characterized by the existence of hubs, i.e. nodes with 
particularly high degrees which act as concentrators of 
connections. Such a finding rules out the possibility that 
the channel system could be organized as a regular net- 
work or a tree, both characterized by similar node de- 
grees throughout. The conditional node density of the 
bone channel structure is shown in Figure |2Ib) . It is 
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clear from such a density that nodes with degree 2 tend 
to connect with nodes with degree 4, which may be re- 
lated to the fact that the circulation of the arteries inside 
the channels follows an Eulerian network, which is char- 
acterized by nodes with even degree |12i] . 

In order to further understand the topological struc- 
ture of the channel system, we considered the recently 
introduced concept of hierarchical clustering coefficient. 
Take a generic node in the network, represented by i. Let 
H{i, d) be the set of nodes which are at distance d from 
node i (the distance between two nodes i and j corre- 
sponds to the number of edges along the shortest path 
between those nodes) . Observe that the consideration of 
subsequent values of d defines a hierarchy around the ref- 
erence node i. Now, the hierarchical clustering coefRcient 
of node i at distance d can be defined as follows 



HCC{i,d) 
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where E{S) is the number of edges between the nodes 
in a set S and \ \S\ \ is the number of elements (or cardi- 
nality) in that set. Observe that the traditional cluster- 
ing coefRcient is obtained by making d = 1 in the above 
equation. It can be verified that the hierarchical clus- 
tering coefficient of node i at distance d quantifies the 
connectivity between those network nodes which are at 
distance d from i. Observe also that < HCC{i, d) < 1, 
with the null value indicating total absence of connec- 
tions between the nodes in H{i,d), while the unit value 
indicates that each node in H{i,d) is connected to all 
other nodes in that set, except itself (loops are not taken 
into account in this work) . Of particular interest is the 
fact that higher clustering coefRcient at a specific dis- 
tance d implies the presence of many cycles composed of 
2d+ 1 edges. 

Figure |2Ic) shows the average value {HCC) (consider- 
ing each network node) of the hierarchical clustering co- 
efficient for d = 1 to 15. Starting at 0.008 for d = 1, the 
values of {HCC) then define a plateau extending from 
d = 2 to 5 and then falls steadily. Such a result provides 
important insights into the topological organization of 
the channel system. First, the low clustering coefficient 
observed for d = 1 expresses the fact that the neighbors 
of each node are weakly connected one another. This in- 
dicates that the channel structure does not involve cycles 
at such a scale. At the same time, the relatively higher 
values of hierarchical clustering coefficients observed from 
d = 2 to 5 implies the existence of several cycles of 5 to 



11 edges and provides further indication that the anlysed 
channel system is not a tree. 

Generally, the arterial vascularization of the major- 
ity of tissues and organs takes place through dichotomic 
bifurcation, i.e. an artery bifurcates into two arteries, 
which in turn bifurcate into four arteries, and so on. The 
obtained experimental results provide cogent indication 
that the vascularization of the bone structure does not 
follow such a general rule. A possible explanation for the 
deviation from the traditional tree and the presence of cy- 
cles could be related to different constraints imposed by 
specific tissues over structural reconstruction after even- 
tual artery obstruction, which can deprive whole regions 
of proper blood flow. While in tissues other than bones 
the organism is known to initiate a neovascularization 
intended to provide bypasses to region irrigation, such 
schemes are virtually impossible in the short term in 
bones, which would demand a complex and long reorgani- 
zation of both Havers and Volkmann channels containing 
the arteries (recall that the bone matrix is rigid) . There- 
fore, it seems that the particularly intense connections 
between nodes from the second to the fifth hierarchical 
levels, as indicated by the obtained hierarchical cluster- 
ing coefficients, represents a possible redistribution sys- 
tem intended to immediately cope with eventual artery 
obstruction. 

All in all, this work has described how concepts 
from three-dimensional reconstruction and image anal- 
ysis, used jointly with state-of-the-art complex network 
theory, can be effectively applied in order to quantify 
certain topological features of biological interconnecting 
systems, particularly the channel distribution network in 
mammals bones. We have verified that such structures 
involve the presence of distribution hubs, present a ten- 
dency for node connections with degree 2 and 4, and 
contain cycles with typical lengths ranging from 5 to 11 
edges. Important biological implications of such topolog- 
ical features have been identified and discussed. Further 
related investigations include the possibility to quantify 
the width and length of each channel in order to infer ad- 
ditional insights about the bone channel system by using 
weighted and geographic complex networks. 
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FIG. 2: The dilog plot of the node degree density, shown in (a), suggests a scale-free network model and indicates the presence 
of distribution hubs. The conditional node degree density is presented in (b), expressing the predominance of connections 
between nodes with degrees 2 and 4. The average values of the hierarchical clustering coefficients, shown in (c), is characterized 
by a plateu extending from d = 2 to 5. 
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